Stability of Weyl metals under impurity scattering 
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We investigate the effects of bulk impurities on the electronic spectrum of Weyl semimetals, a 
recently identified class of Dirac-type materials. Using a T-matrix approach, we study resonant 
scattering due to a localized impurity in tight binding versions of the continuum models recently 
discussed by Burkov, Hook, and Balents, describing perturbed four-component Dirac fermions in 
the vicinity of a critical point. The impurity potential is described by a strength g as well as a 
matrix structure A. Unlike the case in d-wave superconductors, where a zero energy resonance can 
always be induced by varying the impurity scalar and/or magnetic impurity strength, we find that 
for certain types of impurity (A), the Weyl node is protected, and that a scalar impurity will induce 
an intragap resonance over a wide range of scattering stength. A general framework is developed to 
address this question, as well as to determine the dependence of resonance energy on the impurity 
strength. 



Pathbreaking discoveries in the areas of graphene and 
topological insulators have focused attention on a new 
class of materials, known as Dirac materials The 
hallmark of these systems is the existence of one or 
more symmetry-protected Dirac nodes in the electronic 
band structure, where the density of states (DOS) be- 
comes vanishingly small, and which have topological im- 
plications for the existence of gapless chiral edge states 
[4]. Such gapless bulk (three-dimensional) band struc- 
tures have been considered in a number of recent in- 
vestigations [5HTT]. A more general taxonomy was ad- 
vanced by Burkov, Hook, and Balents (BHB) [12 , who 
described the effects of various homogeneous perturba- 
tions on a 3 + 1-dimensional system with a massive 
Dirac point described by the four-component Hamilto- 
nian Hq = Yl't=i + mT^ written in terms of Dirac 
matrices. This provides a minimal model of a system 
with both time-reversal (T) and inversion (X) symme- 
tries, and with a single tuning parameter m which dis- 
tinguishes the massless Dirac point from the normal and 
topological insulating phases when m ^ 0. BHB found 
that for sufficiently strong T or X breaking perturbations, 
an intermediate Weyl semimetal phase generally arises, 
in which the electronic structure is gapless and charac- 
terized by point or line nodes. 

In this paper we consider the effect of localized impu- 
rities on the electronic structure in the Weyl semimetal 
(WS) phase. In particular, we are interested in the occur- 
rence of resonances in the vicinity of zero energy, where 
the density of states vanishes as cj^ in the WS. We say 
that the energy node of a Dirac/ Weyl material model 
H^^"^ is stable if it does not result in a low energy reso- 
nance in the presence of a local impurity potential gAS{x) 
with arbitrary where g is coupling strength and A is 
a matrix describing the scattering potential. Such reso- 
nances give rise to sharp peaks in the density of states 
(DOS) which disrupt the pristine Dirac spectrum [13j . 



In addition to its characteristic low energy electronic 
structure, the WS is also characterized by its topolog- 
ical properties, which interpolate between those of the 
non-topological insulator (NI), where the bulk and edge 
spectra are both gapped, and the topological insulator 
(TI) or - when T is broken - the Chern insulator (CI), 
where the bulk is gapped but gapless edge states par- 
ticipate in quantized surface transport. Bulk transport 
consequences for scalar impurities were considered in refs. 

[HIE!. 

We will focus on the stability of a 4-band tight binding 
model of Weyl materials under local impurity scattering. 
One might guess, based on the more familiar single-band 
problems, that an impurity resonance or bound state can 
be induced at arbitary energy, i.e., no energy is stable. 
This is not true for Weyl materials. Instead, we find that 
stability depends crucially on the type of impurity, which 
mathematically can be classified by its commutation re- 
lations with the F matrices appearing in the local Green's 
function. Typically an impurity is a foreign atom or lo- 
cal crystalline defect in an otherwise pristine material. 
Thus, the impurity potential should always involve a lo- 
cal scalar scattering component. We find that potential 
scattering (local chemical potential on-site change) will 
induce an intragap resonance and therefore will break 
stability at a single particle level. We will present a gen- 
eral framework to address the existence of impurity res- 
onances and bound states, and the dependence of their 
energies on impurity strength. We will illustrate this first 
with a simpler case where both time reversal (T) and 
inversion (X) symmetries are conserved in the impurity- 
free system. There, the (fine-tuned) Dirac node is found 
to be unstable with X-even impurities, but stable with 
X-odd ones. The same approach can be used when T 
and/or X are broken by a homogeneous term rjT^^ ^ and 
the Dirac node is replaced by a pair of Weyl nodes (or a 
line node) for sufficiently strong r]. In these cases, stabil- 
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ity depends not only on the impurity type, but also on r^, 
the strength of the symmetry breaking term. Results will 
be presented for the physically motivated Burkov-Balents 
model [7] of alternating topological/normal insulator lay- 
ers in an external magnetic field {r]) along the stacking 
direction. The impurity classification can be found in 
Table |Tj and the stability phase diagram in Fig. |3] The 
critical field strength r]c is found to be related a kind of 
band inversion. 

Weyl material model - The models we study are lattice 
versions of the continuum models discussed by BHB [12], 
and are defined by the following A;-space Hamiltonian in 
the F-matrix basis, 

3 

H^^\k)=^{k)l^^d,{k)r^m{k)T^^r]r^'' (1) 

i=l 

where ^(/c) = —2t'^- cos ki — Sq, di{k) = — 2tisin/c^, 
m{k) = — 4t'X]i(l ~ coski) — A, and r] is taken to be 
/c-independent for simplicity. We adopt the following F 
matrix convention: F* = (8)<j* {i = 1, 2, 3), F^ = (8)1, 
F^ = -ry 1, and F^^ = z[F^,F^]/2. r' and a' are 
two sets of Pauli matrices acting on the orbital and spin 
degrees of freedom, respectively. In the model, t and 
f are hoppings between same orbitals, and ti is the 
(spin- mixing) hopping between different orbitals. For 
A = 7^ = 0, the model is finely tuned, and there is a 
confluence of the four bands at the single Dirac node 
k = with energy E = —6t — Sq. Nonzero A opens up a 
gap. With 7? = 0, the model is both time-reversal and in- 
version symmetric. Time reversal is defined as T = JCIZ^ 
where JC is complex conjugation and IZ = I <^ io-y . In- 
version is defined as X = F^ (the two "orbitals" being 
opposite inversion eigenstates) . The model will exhibit a 
semimetal phase over a range of r] values. 

T-matrix and single impurity - The effect of localized 
impurities can be studied in the standard T-matrix for- 
malism J5l |16] . We briefly recall the procedure here to 
establish notation. The Green's function of the Hamil- 
tonian H = i^^o) + y is G = G(o) + G^^^TG^^^ where 
G^^\z) = {z- and T = V{I- G^^^V)'^ is the 

T-matrix. Assume H^^"^ is translationally invariant, and 
the impurity potential is localized at the spatial point 
r = 0: Vrr' = g^^r,o^r' ,0- Then the Green's function 
connecting r and r' is Grr' = + G^o^TqqGq^-^,, where 
^00 (^) = (^"^A"^ - g[)o^(^))-^ is the local T matrix, 
and Gqq{z) = ~ ^k"^)'^ unperturbed lo- 

cal Green's function. Here H^^^ is the Fourier transform 
of H^^"^ and N is the number of k points. 

The impurity effect at energy uj can numerically be 
characterized by the matrix norm ||Too|| = \/'^a l^^p 
where A^ are eigenvalues of Too(co' + iO+). For instance, if 
||Too(co'+iO+)|| is large, the local density of states (LDOS) 
p{uj) = —ImJr Grr{^ + ^0^) will in general deviate sig- 
nificantly from the unperturbed one close to the impu- 
rity, yielding a sharp resonance peak. In Fig. [T] we plot 
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FIG. 1: (Color online) Background image: log ||Too(co' + ze)|| 
for T and X symmetric case (77 = 0). Darker color corresponds 
to stronger impurity effect. Unperturbed DOS is shown at 
the bottom. Colored lines interpolating the dark curves are 
obtained by replacing G'^q with ^qq'* (see text). For A = 
F^ (a), the Dirac node is stable. Panel (b) shows results 
for A = F^^ (all curves) and for A = I (blue dashed curves 
only); the Dirac node is unstable. Parameters used are t — 
0.05, ti = -0.5, = -0.25, A = 0,so = -0.3, and lattice size 
Nx — Ny — Nz — 50. Spectral broadening e is set to 0.05. 
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FIG. 2: (Color online) Resonance LDOS at nearest neigh- 
bor site [x^y^z] = [0,0, 1] by scalar and magnetic impurities. 
Solid curves: scalar impurity (A = I). Shaded curves: mag- 
netic scattering (A = F^^). LDOS with different g couplings 
are plotted on different base lines. Gray dotted curve: unper- 
turbed LDOS rescaled. Same parameters as used in Fig. ^ 



log 1 1 Too 1 1 for three impurity forms, showing two qual- 
itatively different cases: in (a) the Dirac node is stable 
because Too has no poles near = 0. For scalar or purely 
magnetic impurities (b), the stability of the Dirac node is 
disrupted by low energy resonances. Fig. [2] plots for case 
(b) the LDOS on the nearest neighbor site of the impu- 
rity. The positions of the resonance peaks agree with the 
large log | |Too 1 1 line in Fig. [ijb). 

Bound states and resonances - Bound states appear 
as poles in G'(cj), which is to say zeros of Tqq{(jj) = 
^-1^-1 _ Q^^^(^i^y This may occur when uj lies outside 
the bulk bands. Resonances are solutions where uo is 



3 



A 


commutation with 
p4 pl2 p4pl2 


g-' yielding det{A/g - 0^°^) = 


physical types 


I p4 pl2 p35 
pl3 pl5 p23 p25 


(+,+,+) 
(+,-,-) 
(-,+,-) 
(-,-,+) 


Sy/g = a + S4&1 + 51262 + 563 
= (a + 5461)^ - (62 + 5463)^ 
1/3^ = (a + 51262)^ - (61 + 51263)^ 
= ia + sb3f -{bi+sb2f 


scalar, Jz and/or layer-swapping 
J|| or with layer- swapping 
layer- mixing or with Jz 
J|| with layer-mixing 



TABLE I: Impurity classification for point-node Weyl semi-metal, with F^^ = F^^ in eqn.[l] In the second column, + denotes 
commute and — anticommute. a and bi are defined after eqn. |3] principal values are implicitly taken. S4,si2 and sa are 
eigenvalues of F'^,F^^ and A, respectively, and take the value ±1. s = S4S12. See Fig. [sjfor stability phase diagram of the last 
three classes. 



extended to have a finite imaginary part, which is in- 
versely proportional to the resonance lifetime. In solving 
numerically for the resonances, it is convenient to de- 

fine g^l\uj) = H^oo + + [G^o + ^e)]^), where 
e ^ is always greater than the energy spacing between 
consecutive quantized levels. We then seek solutions to 
Ua{uj) = g~^^ where {ua{uj)} are the eigenvalues of Gqq 
This prescription works well when the bulk DOS is small, 
as we show in Fig. 1, and is exact when the DOS van- 
ishes. Thus, if, for a given real value of cj, at least one of 
the eigenvalues {ua{oj)} is real, and if the bulk DOS is 
small, then a bound state or resonance exists for coupling 
g = l/ua{oj). If all Ua{oj) have imaginary part, then uj 
is a stable energy. An immediate consequence is that an 
impurity with [A,^qq^] = will for some g disrupt the 
stability at cj, because product of commuting Hermitian 
matrices has real eigenvalues. Single-band problems fall 
in this category as A = I. 

Dirac node - To illustrate the stability criteria, we 
consider first the case with ?7 = 0, where both X and T 
symmetries are present. Inverting Eq. [l] one finds that 
the only F matrix appearing in G^J is F^, which is also 
the inversion operator, 

with a = ^ E - m]/D{k), ^ = ^ Efc m{k)/D{k), 
and D{k) = [uj- ^{k)]^ - ^. d^^{k) - m^{k). In the fol- 
lowing, we will restrict to impurities where A is either the 
identity or a single F-matrix; linear combinations thereof 
can be analyzed in the same fashion. A belongs to one 
of two classes: F^ which commutes with F^, and F^ 
which anticommutes. Following our general criteria, we 
solve for real eigenvalues of ^qq^A. The commuting case 
has already been discussed, thus inversion-even impuri- 
ties may disrupt stability at any energy, including the 
Dirac node. As an example, the purely magnetic impu- 
rity A = F^^ = I(g)<j^ is shown in Fig.jljb), where dashed 
lines are the eigenvalues Ua- The corresponding LDOS 
on the nearest neighbor site to the impurity is plotted 
in Fig. [2] for various g which induces resonance around 
the Dirac node. For inversion-odd impurities, we have 
a^o^F^ = aF^ + ^F^F^, and u{uj) = ±^a{ujY -h{ujY, 



where principal values of a and h are implicitly taken, 
as will be all such coefficients in the rest of the paper. 
Reality oi g~^ = u then requires \ci{uo)\ > \b{uj)\. A rep- 
resentative case with A = F^ = is shown in 
Fig. [TJb), where real u are the dashed blue lines. 

The region between the two LDOS "towers" is stable 
in Fig. [ijb). This in fact is true for all inversion-odd 
impurities and is not accidental. This follows analytically 
from the band center approximation (BCA) in which the 
local Green's function, ^qq^ , is replaced with the Green's 
function of the local Hamiltonian, G{uj) = {uj — Hqq)~^^ 
where Hqq = -^^j^H^^\k). It is easy to verify that 

Hqq has eigenvalues Q± = —Eq ± with 6 = \12t^ + A|, 
which are doubly degenerate due to inversion symmetry. 
We shall call Q± band centers as they represent in some 
sense the average position of the bands. The eigenvalues 
of G{uj) are thus l/{uj — ^±). Comparing with eqn. [2j 
we find the BCA for a and b as a^^^ = ^^_^^'^2_§2 and 
j^BCA __ Invoking our earlier results, g — 

is imaginary for uj G [1^-, Thus the Dirac node, 
being in between the band centers in general, is stable 
with X-odd impurities. 

Weyl semimetal phase - BHB [12] found a WS phase 
intervenes when T or X is broken sufficiently strongly. 
Depending on F^^ in Eq. [l] the two central bands may 
touch zone either at discrete points or along a line in the 
Brillouin zone. We now consider the effect of local im- 
purity on the WS phase induced by the T-breaking term 
F^^ = F^^ = I (g) (7^. A complete classification will be 
presented elsewhere [17 . The Hamiltonian Eq. [ijcan be 
used to model an alternating stack of topological insula- 
tor and normal insulator layers immersed in a magnetic 
field along the stacking direction (the F^^ term), as pro- 
posed in Ref. i The local Green's function is 

G^^Jiuj) = a{uj) I + bi{uj) F^ + b2{uj) Fi2 + i,^^^) 

(3) 

where all coefficients can be obtained analytically fT7]. 
Note that all four terms in Eq. [3] mutually commute. 

Following the general stability criteria, we find that im- 
purities fall into four classes, based on the commutation 
relations of A with the three F matrices in eqn. [3j The 
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FIG. 3: (Color online) Stability of T-breaking Weyl 
semimetal. The vertical axis 77 is the strength of the 
T-breaking term. Red: stable zone for impurity classes 
(+,—,—) and (—,—,+). Blue: stable zone for impurity 
classes ( — ,+,—) and ( — ,—,+). See Table [l] for impurity 
classification. Solid black lines mark the two band edges 
bounding the central gap. They touch from around 77 = 1 
to 5, corresponding to the Weyl semimetal phase. The black 
lines are broken around 77 = 1 due to the closing of the in- 
direct gap. Dotted gray lines are the stable zone boundaries 
given by the band center approximation. Parameters used are 
t = 0.05, ti = -0.5, = -0.25, A = l,so = -0.3, on a lattice 
of = Ny = = 50. 



results are summarized in the Table [Tj The fully commut- 
ing class (-h, +, +) always yields a real g solution for any 
uj and T]. Each of other three classes yields a stable zone 
in the uj-t] plane where g is imaginary. A phase diagram 
is shown in Fig. [3j where red is the stable zone of the 
class (+, — , — ), blue of (— , +, — ), and (— , — , +) is stable 
in both red and blue zones. The general shape of the 
zone boundary can be understood by invoking the band 
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center approximation again, viz. Gqq {(^) — {00 — 
This yields the two dotted lines in Fig. [3] given by 

uj = ^± = -eoT{\5\-\n\) , (4) 

which are nothing but the central two of the band cen- 
ters (eigenvalues of Hqq). A band center inversion oc- 
curs when = at which point Q-^ = It may 
be viewed as a stability critical point (critical external 
magnetic field) in the sense that the stable zones of the 
classes (+,—,—) and (—,+,—) are each restricted to only 
one side of rjc One may think of impurity levels as form- 
ing an energy band parameterized by the strength g^ then 
stable zones in the phase diagram are gaps of the "im- 
purity band" , and 77c marks a transition between gapless 
and gapped phases. This is reminiscent of the Bloch band 
inversion associated with the topological phase transition 
in Chern insulators. 

Average T -matrix approximation (ATA) - The reso- 
nance/bound states as discussed before will broaden into 
a "band" when multiple impurities are present. Here we 



3.0 

2.5 

2.0 

O 1.5 
Q 

1.0 
0.5 
0.0 



0123456789 

i\ 

.^Jl^^ 


, clean 

i ATA 


""""^ ./I 1 

/ \ / M 
/ V '^J 





I " " I " " I " " I " " I " " I " " I " " I " " I 
-10 -8 -6 -4 -2 2 4 6 8 10 



1.2 
0.6 

-0.6 
-1.2 



FIG. 4: (Color online) Average DOS using ATA for F^^ = F^^ 
(external magnetic field along z) and scalar impurity A = I. 
Impurity concentration is c = 10% and impurity strengths are 
uniformly distributed in ^ G (0, 10]. Solid blue line: ATA re- 
sult. Dash-dot black line: DOS of clean system. Background 
color image: log||Too|| of single impurity as a function of uj 
and , similar to Fig.hJ Dotted horizontal line marks maxi- 
mum g used in ATA (^ = 0.1), the ATA DOS is significantly 
enhanced in the range of u in which the high log ||Too|| lines 
exist above g~^ = 0.1. 77 is set to 3, putting the clean sys- 
tem in the WS phase (Fig.|3|. Other parameters used are the 
same as in Fig. |3] 



consider an ensemble of local impurities of the same ma- 
trix form, with a homogeneous spatial concentration c 
and a distribution of strength f{g). We find that av- 
erage DOS is enhanced at resonance energy a;res(^) ob- 
tained from a single impurity when f{g) is relatively high. 
In the ATA formalism [151 E] 5 translational invariance 
is restored after statistical averaging, and the effect of 
the impurity ensemble is captured by the local self en- 
ergy Sloe = c(Too)[l + cG^'o^(Too)]-\ where (• • • ) de- 
notes averaging over f{g). The local Green's function 

Eioc) and the aver- 
z0+). In Fig.ji) we 



is Gioc(^) = ^Efcl/(^-^'°nfc) 
age DOS is p{uj) = — lmTrGioc(^ 

show the average DOS of a sample immersed in an exter- 
nal magnetic field (F^^ = F^^) with magnetic impurities 
breaking local inversion symmetry (A = F^). The im- 
purity strengths are uniformly distributed in ^ G (0, 10]. 
The difference in DOS with the clean system agrees with 
the single impurity | |Too| | (shown as colored background). 

Conclusion. We have investigated the stability of the 
electronic spectrum of the Weyl semimetal phase with 
respect to local impurity scattering. 
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